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Abstract. Hanle-effect observations of forbidden coronal line transitions and recently also 
longitudinal Zeeman-effect measurements of coronal lines show quantitative signatures of 
the weak coronal magnetic field. The interpretation of these observations is, however, com- 
plicated by the fact that they are the result of line-of-sight integrations through the optically 
thin corona. We study by means of simulated observations the possibility of applying tomo- 
graphic techniques in order to reconstruct the 3D magnetic field configuration in the solar 
corona from these observations. The reconstruction problem relates to a family of similar 
problems termed vector tomography. It is shown that Zeeman data and Hanle data alone ob- 
tained from vantage points in the ecliptic plane alone are sensitive only to certain magnetic 
field structures. For a full reconstruction it is necessary to combine the longitudinal Zeeman 
and Hanle effect data. 
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1. Introduction 

The coronal magnetic field is the main driv- 
ing force for most plasma processes in the in- 
ner corona. To understand the physics of the 
corona a detailed knowledge of the state of 
the coronal field is therefore essential. Unlike 
the photosphere, however, the low density and 
high temperature of the coronal plasma make 
direct field measurements in the corona to be 
rather difficult. Conventionally, the magnetic 
field of the corona is therefore estimated by 
means of extrapolations from its photospheric 
boundary values. Until now, potential field ap- 
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proximations as an extrapolation model were 
quite popular. The omission of all electrical 
currents in the corona, however, misses an im- 
portant part of coronal magnetic field driven 
physics: as the potential field is the lowest en- 
ergy state of the field wi th respect to normal 
boundary conditions (e.g. lSakuraill989f) it can- 
not account for the energy stored in the field 
which is released in dynamical processes such 
as flares and CMEs. 

With the advent of vector magnetograms 
from the solar surface it became possible to 
employ more a sophisticated and more realis- 
tic model assumption, the force-free condition, 
for the field extrapolation. There is currently 
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a considerable interest in algorithms which 
solve this non-linear and ill-posed extrapola- 
tion problem (see, e.g., the contribution by S. 
Regnier in this issue). It is obvious, however, 
that extrapolations can only reproduce those 
structures of the magnetic field which have a 
measure-able imprint on the surface field. The 
results are deemed to become less reliable the 
larger the distance from the surface. 

In this paper we investigate therefore 
whether spectropolarimetric coronal measure- 
ments of the magnetic field can be used for 
its reconstructio n. These meas urements have a 
long tradition. ICharvinl (|T965) suggested more 
than 40 years ago to measure the Hanle effect 
in the corona. One of the first to perform ac- 
cording observa tions of the cor onal green line 
at 530.3 nm was lArnaudl(ll982l) . But full spec- 
tropolarimetric measurements of coronal lines 
are still today a challenge. For a long time, 
the observation of the longitudinal Z eeman ef- 
fect was out of reach. Only recently [Lin et alJ 
(2000,|2004|) could demonstrate that fuU spec- 
tropolarmetric observations in the corona can 
be achieved. This success is partly due to im- 
proved instrumentation but also profits largely 
from a change to infrared lines where the 
Zeeman split is increased relative to the ther- 
mal line broadening. The Zeeman observations 
of the Fe XIII 1075 nm line, however, still re- 
quired almost an hour of integration time. 

For forbidden coronal lines with a life time 
much larger than the Larmor period of the ex- 
cited atom a first order interpretation of the ob- 
served spectropolarimetric signals can be given 
as follows. The Stokes V signal is proportional 
to the line-of-sight magnetic field component 
B]l at the location where the line is emitted, 
the Stokes Q and U polarisation signals indi- 
cate the orientation B±/B ± normal to the line 
of sight. 

This picture of the line formation, how- 
ever, is a strong simplification. The physics 
of the emission of these forbidden lines has 
been studied in great detail in recent years , 



e.g bvlHousd dl977h . ISahal-Brechol d 19771) 
and|0uerfeldj(l982) and our understanding has 
been largely refined with respect to the above 
crude interpretation. 



The above crude picture, however, also a 
neglects the fact that the coronal observations 
are line-of-sight integrals through an optically 
thin medium and do not represent directly its 
local properties. The goal of our study is to ob- 
tain a more sophisticated interpretation of the 
observations also in this respect. In particular, 
we want to investigate the whether these ob- 
servations suffice to determine a global model 
of the coronal magnetic field. Due to a lack 
of space we can here only briefly discuss the 
approach we have adopted (chapter |2| and 
present initial results of some of our test cal- 
culations (chapter 0. Necessary future work 
work is outlined in the final s ection (chapter 
0). More details can be fo und in Kj-ama] (I2005J) 
and lKramar et all (120061) . 



2. The inversion problem 

The inversion of line-of-sight observations is 
commonly referred to as tomography. In many 
fields this is a well established technique which 
solves for the distribution of an isotropic scalar 
quantity in a limited region of space from line- 
of-sight integrals of the scalar along a dis- 
crete number of directions. The inversion is ill- 
posed but has no null space which means that 
all structures of the density distribution can in 
principle retrieved, the quality of the recon- 
struction depends on the resolution and signal 
to noise ratio of the individual images taken 
and on the angular sampling of the image view 
directions. 

Tomog raphy has al so been applied to solar 
physics. Davila i 19941) proposed to use coro- 
nagraph observations to construct the coro- 
nal plasma density. A ccording calcu lations 
have been performed by Izidowitj i 19991) and 
iFrazin and Janzenl d2002l) . For the application 
of tomography to space observations some 
additional difficulties have to be coped with 
which are absent in laboratory applications. 
The observations are usually made with a coro- 
nagraph with the Sun being occulted. This 
causes a data loss at the center of each image 
which increases the condition of the inversion 
problem. Observations are traditionally made 
from ground-based telescopes or from space 
craft close to the ecliptic more or less contin- 
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Fig. 1. Illustration of the inability to see divB and e.ios xcurlB from certain line of sight integrals 



uously over half a rotation period of the Sun. 
The tilt of the Sun's rotation axis with respect 
to the ecliptic and the instationarity of coronal 
structures during the time of the observation 
may seriously affect the quality of the recon- 
struction. 

The inversion problem we are faced with 
bears a further very fundamental complication. 
We want to reconstruct a field rather than a 
scalar quantity and as observations we have 
line of sight integrals of the field components 
(or functions thereof) along and perpendicular 
to the view direction. This extension to scalar 
tomograph y has been in vestigated in recent 
years (e.g. lHowardf T996) and has been termed 
vector tomography. If we adopt for a moment 
the very simplified interpretation of the spec- 
tropolarimetric observations mentioned in the 
previous chapter, for the Hanle effect we for 
the moment even assume that we can measure 
the strength of the field in the plane of the sky, 
we would observe quantities like 



LOS LOS 



d€ (1) 



It is a well known result from vector tomogra- 
phy, that the top integral is completely insen- 
sitive to the any divergent part of the field B 
and the bottom integral does not respond to the 
components of the curl of B normal to the line 
of sight (see Fig. [TJ. Hence we have to deal 
with an extended null spaces in the field we 
want to reconstruct from the above measure- 
ments. 

How far the problem of a finite null space 
also applies to real Zeeman and Hanle ef- 
fect observations is not obvious because their 
emissivities are more complicated than was as- 
sumed in (H}. We therefore have developed 



an inversion procedure to numerically test 
the condition of the field inversion. In our 
test comp utations we use the emissivity ex- 
pressions (lHouselll977t ISahal-Brecho3ll977t 
lOuerfeldll 19821) 
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Here, £ and A are proportional to the ion den- 
sity and to the sum and difference, respectively, 
of the upper sublevel population of the emit- 
ting electron transition, g is the effective Lande 
factor of the transition. The wavelength depen- 
dence of the different Stokes lines is neglected 
here assuming that representative moments of 
the line signal were taken. The magnetic field 
enters into (|2} through the field intensity in the 
Larmor period o>l and the field angles with 
respect to the line of sight and a, its orienta- 
tion normal to the line of sight. 

Our inversion procedure is based on a least 
square iteration between the polarimetric ob- 
servations and a forward modelling of the data 
calculated from the respective line of sight in- 
tegrals of (0 assuming a model field B. The 
field is then successively improved until the 
difference reaches the instrument's noise level. 
According to the above arguments we have to 
expect that the forward modelling may have 
null spaces and our minimisation may not yield 
unambiguous results. To stabilise our inver- 
sion, we therefore add a second term to our ex- 
pression L(B) to be minimised 
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+ J |divB| 2 dv (3) 

corona 

This second term is an obvious constraint to 
the magnetic field and it has far reaching con- 
sequences. The calculation of the divergence at 
the inner coronal boundary requires the normal 
component of B on the Sun's surface as bound- 
ary condition. Hence our data set has to be ex- 
tended to include not only the coronal observa- 
tions but also the photospheric surface magne- 
tograms. 

We consider this extension an advantage 
because this additional information is merged 
in a natural way with our inversion procedure. 
If, e.g., we had no coronal observations at all, 
the procedure would only minimise the diver- 
gence term in d3j. If in addition we could in- 
sure that the total field energy J \B\ 2 dv is also 
minimised, the solution will be the potential 
field which complies with the measured sur- 
face magnetogram. The energy minimisation, 
however, is implicitly taken care of if a proper 
minimisation algorithm (e.g., conjugent gradi- 
ents) is chosen along with a minimum energy 
initial field for the iteration. Hence we obtain 
the potential field approximation of the coro- 
nal magnetic field for free and every individual 
coronal observation included in © drives our 
solution towards a more realistic, current car- 
rying corona. 

The problem we have to investigate can 
be restated much more precisely now: is the 
first term in ©, i.e. the coronal spectropo- 
larimetric data, sufficient to resolve the mag- 
netic field perturbations of all current systems 
to be expected in the corona ? Reversely, are 
there current systems which the coronal obser- 
vations do not respond to and which will con- 
sequently not appear in our solution. We note 
in passing that we could easily constrain the 
magnetic field even further by adding a third 
term in (0 proportional to the integral of the 
squared j x B force in the corona. This would 
additionally stabilise our inversion and it is in 
fact exactly the approach some popular surface 
field extrapolation schemes make use of (see 
S. Regnier's contribution in this issue), except 
that there the term of the coronal observations 
in 0]) is not accounted for. 



3. Test of the inversion approach 

In this section we report on first initial test cal- 
culations using the approach described above. 
The test consisted of the retrieval of given 
coronal magnetic field configurations from 
simulated observations. The coronal density 
was assumed to fall o ff radially acco rding to 
a classical power law dNewkirkll 19701) and as- 
sumed was known for the field inversion be- 
cause it can, at least approximately, be ob- 
tained from a scalar tomography inversion of 
the line intensity. 

In Fig.|2]we show the first model we tested. 
It consists of a dipole field with an isolated 
current loop across the equator confined on a 
meridional plane. The simulated data set was 
generated form this model field as if observed 
by a space craft on an ecliptic plane tilted by 
10 degrees with respect to the z axis. The data 
set thus comprised images of the Stokes signals 
taken from 36 equidistantly distributed view- 
ing directions with 5% noise added to the cal- 
culated signals intensities. 

The inversion cannot avoid to reconstruct 
the field in the whole domain of the corona. 
This should not be a problem where the field 
is a potential field and we therefore empha- 
sise the perturbed field region when present- 
ing the inversion results. Fig. [3] shows these 
field perturbations in the equatorial plane in the 
immediate neighbourhood of the point where 
the current loop intersects the plane. Note that 
the background dipole field is normal to this 
plane so that it does not show in Fig. [3] For 
test purposes we have here only made use of 
either the Zeeman (centre panel in Fig. [3j or 
the Hanle (right panel) observations for the 
reconstruction. For comparison, the left panel 
displays the perturbation of the original field. 
Obviously, the Zeeman observations yield a 
much better reconstruction than the Hanle data. 
The reason seems to be that for this orientation 
of the current, the Hanle data responds very lit- 
tle while the Zeeman effect alone provides suf- 
ficient information to resolve this kind of field 
perturbation. 

As an alternative model, we investigated a 
dipole field with a current loop in the equato- 
rial plane. Again in Fig. [3] we show the origi- 
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Fig. 3. Field perturbations in the equatorial plane close to where the current loop of model in Fig.[2]inter- 
sects the plane. From left to right we show the original field, the reconstruction taking only account of the 
Zeeman observations and of only the Hanle observations in l[3j. 




Fig. 4. Coronal magnetic field model 2 with a current loop in the equatorial x, y plane. We also show the 
distorted field lines. 



nal and the reconstructed magnetic field in the 
x, y-plane normal to the current loop and con- 
centrated to the neighbourhood of the current 
intersection. Note that now the background 
dipole field is superposed and the effect of the 
current cannot so clearly be discerned. It can 
be concluded though that now the Zeeman ob- 
servations miss the field perturbation entirely 
and only allow to reconstruct the background 



dipole field. The Hanle observations on the 
other hand give a decent response and its in- 
formation included in the inversion reproduces 
to some extent the effect of the current. 

4. Summary and outlook 

We have presented first results which aim to 
use coronal spectropolarimetric observation to 
reconstruct the magnetic field of the corona. 
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Fig. 5. Field perturbations in the meridional plane centred between the loop's foot points and close to where 
the current loop of model in Fig.[2]intersects the plane. From left to right we show as in Fig.[3]the original 
field, the reconstruction taking only account of the Zeeman observations and of only the Hanle observations 
in (3). 



From our findings we can state that Hanle or 
Zeeman observation alone are not sufficient for 
the reconstruction. Instead the Stokes V, Q and 
U components and in addition surface magne- 
tograms are necessarily required to be sensitive 
to the most obvious coronal current systems. 

We are confident that this data set is also 
sufficient to yield a realistic coronal magnetic 
field model. This, however, has to be verified 
in future experiments. If it turns out that the 
inversion is too unstable to yield satisfactory 
field models, we still have the option, as men- 
tioned above, to include a force-free constraint 
in our minimising function 0]). 

A future critical test is the dependency of 
the inversion output on the noise level of the 
observations. This test besides its insight into 
the condition of the inversion problem can also 
provide us with an estimate of the maximum 
tolerable signal to noise ratio. This way we 
may also optimise the tradeoff between image 
noise and angular resolution since the Zeeman 
observations require long integration times. We 
are confident that this information could be 
quite helpful to improve real observations. 

We are just at the beginning with our in- 
vestigations of this new method and there is 
quite some work ahead. In future we intent 
to also include the Stokes I component in our 
data set and treat the particle density as addi- 
tional variable so that an a-priori density model 
is no more required. In the more distant fu- 
ture we may also consider to take account of 



the full spectral profile of the Stokes measure- 
ments and in return relax a-priori assumptions 
about the coronal temperature distrubution. 
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